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The basic objective of this work is to compare and evaluate the performance of 
various turbulence models which are used in the numerical simulation of complex 
turbulent flows. The approach utilizes the Reynolds-averaged compressible Navier- 
Stokes equations in which the Reynolds stresses and heat fluxes are mathematically 
modeled by suitable turbulence models. This paper focuses on the simulation of 
transonic flows about a 2-D airfoil and the estimation of performance character- 
istics such as lift and drag. Attention is given to separated flows where differ- 
ences between model predictions are generally greater than those observed for 
unseparated flows. 

The turbulence models used in the present study are eddy viscosity models which 
include the family of zero-, one-, and two-equation models. Zero- and one-equation 
models are numerically the simplest of the eddy viscosity models but lack the 
generality of the two-equation models because their length scales must be determined 
algebraically rather than from a field equation. However, two-equation models are 
more complicated than zero- and one-equation models, and these complications can 
sometimes lead to numerical difficulties. 

Six turbulence models are studied in this paper. They are the zero-equation 
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models of Cebeci and Smith , Baldwin and Lomax , and Johnson and King , and the two- 
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equation k-e and q-w models respectively of Chien and the present author, Refs. 5 
and 6. These models are listed in Table I. All of these models utilize the proce- 
dure of integration-to-the-wall in which no-slip boundary conditions are applied at 
solid surfaces. For the zero-equation models, a common formulation is used in the 
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airfoil wake region which is due to Cebeci . The location of transition in the 
computations was taken from the experimental boundary-layer trip locations except in 
the cases of the two-equation q-w model predictions where the transition locations 
were allowed to occur naturally. 

The numerical differencing method consists of a second-order implicit upwind 

differencing algorithm combined with a finite volume discretization technique which 
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produce accurate resolution of shock waves and discontinuities . Additional fea- 
tures of the method include the use of spatially varying time steps to speed con- 
vergence to a steady state and an implicit treatment of boundary conditions. 
Inviscid boundary conditions at far field boundaries are based on the method of 
characteristics and account for circulation due to lift. 

Paper presented at 10th International Conference on Numerical Methods in Fluid 
Dynamics, Beijing, China, June 23-27, 1986. 


1 



The airfoil investigated in the present study was the RAE 2822 airfoil which 
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has extensive experimental documentation . Three cases were investigated. Cases 1, 
9, and 10, which correspond respectively to unseparated subcritical flow, unsepa- 
rated supercritical flow and separated supercritical flow. The experimental condi- 
tions of Mach number, Reynolds number and geometric angle of attack for each of the 
three cases are shown in Table I. Since the calculations were done in free air, the 
angle of attack must be changed from the experimental or geometric angle of attack 
to account for wind tunnel wall effects. The angles of attack used in the calcula- 
tions are shown in Table I and were obtained from the recommendations in Ref. 9. 

The numerical grid used in the computations is shown in Fig (1)* It consists 
of an algebraically generated 240x60 C-grid. The mesh spacing in the y direction at 
the surface was such that y + at the first mesh point was less than one everywhere on 
the airfoil. The spacing in x over the central portion of the airfoil was AX/0.013 
with the spacings reduced at the leading and trailing edges. The far field boundary 
was placed approximately 20 chords from the airfoil surface. 

Surface pressure and upper surface skin friction and displacement thickness 
distributions are compared for Case 1 in Fig. (2). It is apparent that the surface 
pressure distributions predicted by the turbulence models are in good agreement with 
one another and with experiment. Larger differences between model predictions and 
experiment are indicated by the skin friction and displacement thickness 
distributions. Models which are in best agreement with experiment on skin friction, 
especially near the trailing edge, are the C-S, B-L and q-w2 models. The best 
models in predicting displacement thickness are the J-K and q-w2 models. The lift 
and drag predictions for this case are shown in table I. It appears that all three 
zero-equation models give better predictions of drag than the two-equation models 
while the J-K and q-w2 models give the best predictions of lift. 

The second case studied was the unseparated supercritical case 9. Computed 
Mach contours for this case obtained using the J-K model are shown in Fig. (3). 
Surface pressure and upper surface skin friction and displacement thickness distri- 
butions are shown in Fig. (4). Computations of surface pressure indicate that all 
models are in close agreement with experiment except near the shock wave where the 
J-K and q-w2 models predict shock locations which are slightly upstream of the 

locations given by experiment and the other models. With regard to skin friction, 

the C-S, B-L and q-w2 models give the best predictions downstream of the shock wave 
while the J-K, q-co2 and k-e models are closest in their predictions of displacement 
thickness. The k-e model shows an unrealistically large increase in skin friction 
downstream of the shock wave. From Table I it is evident that the models giving the 

best overall predictions of lift and drag are the J-K, q-u)2 and k-e models. 

Computed Mach contours and surface distributions for Case 10 are shown in Figs. 
(5) and (6). In this case the flow downstream of the shock wave is separated for 
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some distance. The corresponding model predictions show more differences in surface 
pressure distributions and especially shock wave locations than in the previous 
cases and this gives rise to larger differences in lift and drag. It is evident 
that the J-K model provides the best overall predictions of surface pressure, skin 
friction and displacement thickness distributions as well as lift and drag (Table 
I). It should be noted in this case that the calculations using the J-K and q-o>2 
models were unsteady with the shock wave undergoing small periodic oscillations 
about its mean position. The corresponding oscillations in lift were from 2% (J-K) 
to 4% (q-o)2) of their mean values. The k-e model showed an anomalous behavior (or 
weak instability) in skin friction for this case which is believed to result from 
numerical stiffness associated with the low-Reynolds-number damping terms of the 
model (see Ref. 5). 

Numerical efficiency of the various turbulence models was measured by the 
computing time required to achieve steady state, or in the case of unsteady flows, 
the time required to achieve a periodic state. For the unseparated cases 1 and 9 
about 1000 steps were needed for convergence (to approximately 3 significant Figures 
in C^). Computing time for the zero-equation models was about 480 seconds (on the 
NASA Ames Cray XMP/48) and about 600 seconds for the two-equation models. For case 
10, the models which predicted steady solutions (i.e., C-S, B-L, k-e, q-o>l) required 
approximately the same time as cases 1 and 9. The unsteady solutions (J-K and q-^2 
models) required about 2.5 times more computing time than the other cases and a 
constant (instead of spatially varying) time step was used to maintain time 
accuracy. 

The primary conclusions resulting from this work are that differences between 
experiment and turbulence model predictions which are relatively small for subcriti- 
cal unseparated transonic flows, become greater for supercritical separated flows. 
The principal differences between model prediction of lift and drag for the super- 
critical cases appear to be the result of differences in predicted shock location. 
The best overall model in predicting the three cases was the J-K model. 
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Table I - Turbulence Models and Experimental/Ccmputational 
Characteristics - RAE2 822 Airfoil 


Model 

Originator 

Case 

C L 

1 

C D 

Case 

C L 

9 

C D 

Case 

C L 

10 

C D 

Zero-equation models 







C-S 

Cebeci-Smith 1 

.610 

.0086 

.833 

.0172 

.815 

.0269 

B-L 

2 

Baldwin-Lanax 

.620 

.0087 

.861 

.0185 

.859 

.0298 

J-K 

Johnson-King 3 

.573 

.0087 

.787 

.0159 

.745 

.0243 

Two-equation models 







k- e 

4 

Chien 

.599 

.0094 

.821 

.0179 

.776 

.0268 

q-a) 1 

Coakley^ 

.607 

.0095 

.846 

.0181 

.836 

.0284 

q-w 2 

Coakley*’ 

.576 

.0093 

.783 

.0159 

.726 

.0267 

9 

Experiment 

.566 

.0085 

.803 

.0168 

.743 

.0242 

a(gean) 

, a(comp): 

2.40 , 

1.93 

3.19 , 

2.80 

3.19 , 

2.80 

M co 

, R xlO -6 : 

.676 , 

5.70 

.730 , 

6.50 

.750 , 

6.20 





Fig. 6 Surface pressure, upper surface 
skin friction and displacement thickness 
distributions. Case 10 
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